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5h , 

C^h, We examine the Casimir energy of 5D electromagnetism from the recent standpoint. 

. The bulk geometry is flat. The Z2 symmetry and periodic property, for the extra coordi- 

nate, are taken into account. After confirming the consistency with the past result, we do 
' new things based on a new regularization. In the treatment of the divergences, we introduce 

IR and UV cutoffs and restrict the (4D momentum, extra coordinate)-integral region. The 
regularized configuration is the sphere lattice, in the 4D continuum space, which changes 
along the extra coordinate. The change (renormalization flow) is specified by the minimal 
area principle; hence, this regularization configuration is stringlike. We do the analysis not 
in the Kaluza-Klein expanded form but in a closed form. We do not use any perturbation. 
The formalism is based on the heat-kernel approach using the position/momentum propaga- 
tor. Interesting relations between the heat kernels and the P/M propagators are obtained, 
where we introduce the generalized P/M propagators. A useful expression of the Casimir 
, energy, in terms of the P/M propagator, is obtained. The restricted region approach is 

■ replaced by the weight-function approach in the latter-half description. Its meaning, in re- 
■^j- ' lation to space-time quantization, is argued. Finite Casimir energy is numerically obtained. 
\Q , The compactification-size parameter (periodicity) suffers from the renormalization effect. 

Numerical evaluation is exploited. In particular, the minimal surface lines in the 5D flat 
(f) ■ space are obtained both numerically using the Runge-Kutta method and analytically using 

a general solution. 

OO ' §1. Introduction 

O ' 

As a unified theory of the four forces in nature, higher-dimensional models have 
a long history since the papers by Kaluza 1 ) and Klein. 2 ) The simplest one unifies 
the forces of gravitons, photons and dilatons. The quantu m ef fects are evaluated 
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by Appelquist and Chodos. 3 ** They evaluated Casimir energ v**)| and the result has 
been giving us a standard image of the contraction of the extra space, that is, when 
compactification takes place, the extra space shrinks to a size slightly larger than 
the Planck length. 

Higher-dimensional models have, at present, the defect that they, by themselves, 
are unrenormalizable. In Ref.3), UV divergence appears as quintic divergence A 5 of 
the cosmological term. It simply means that we have no appropriate procedure 
for defining physical quantities within the quantum field theory (QFT). One can, at 
this point, have the standpoint that they are effective theories that should be derived 
from more fundamental models such as the string theory, M-theory, and D-brane. In 
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the present study, we pursue the possibility that there is an appropriate procedure 
for defining physical quantities within the higher-dimensional QFT. We propose a 
procedure and show that it works well. 

There are some approaches to solving the above problem. One of them is the 
deconstruction model. 5 )' 6 ) We discretize the extra coordinate and choose an appro- 
priate finite number of "branes" keeping gauge invariance. This is a commonly used 
approach at present. Some interesting results are reported. 7 )' 8 ) The other approach 
is based on the regularization using a position-dependent cutoff. 9 -* The integral re- 
gion is restricted appropriately. The restriction requirement comes from the analysis 
of propagator behaviour. 9 )' 10 ) Spiritually, the holograpy idea is behind the restric- 
tion procedure. The present motivation comes from the question, "Can we find the 
reason why the restriction process is necessary within the framework of the 5D QFT, 
not using the string theory and related supergravity theories?" 

We introduce a new regularization inspired by the partial success of Randall- 
Schwartz's result. We associate the regularization (in 4D world) cutoffs running 
along the extra axis y, with the minimal area surfaces in the bulk. In this way, 
the stringlike (surface) configuration (closed string) is introduced in the present 
approach. This is quite in contrast to the usual string theory approach. The present 
stringlike configuration appears not from the propagation of strings but from the 
necessity of the restriction of the integral region in the bulk space. 

The original approach to the renormalization flow interpretation of bulk be- 
haviour relies on the AdS/CFT and 5D supergravity. 11 ) -16 ) The present approach 
does not rely on them. We directly use the minimal area principle, the essence of the 
string theory, 17 )~ 19 ) in the regularization procedure. This is new in the development 
of the quantum field theory. 

This paper is organized as follows. In §2, we review the 5D quantum electromag- 
netism in the recent standpoint. Casimir energy is obtained from the KK-expansion 
approach. In §3, the same quantity of §2 is dealt with in the heat-kernel method, 
and Casimir energy is expressed in a closed form in terms of the P/M propagator. 
The closed expression of Casimir energy is numerically evaluated and its equivalence 
with the result of §2 is confirmed in §4. Here we introduce UV and IR regularization 
parameters in the(4D momentum, extra coordinate) space. A new idea about UV 
and IR regularization is presented in §5. The minimal surface principle is introduced. 
The sphere lattice and renormalization are explained. In §6, an improved regulariza- 
tion procedure is presented where a weight function is introduced. Here, again the 
minimal surface principle is taken. The definition of the weight function is given in 
§7. In §8, we present the conclusions. We prepare two appendices to supplement the 
text. Appendix A deals with the analytic solution of the minimal surface curve in the 
5D flat space. Appendix B provides an explanation of the numerical confirmation of 
the (approximate) equality of the minimal surface curve and the dominant path in 
Casimir energy calculation. 
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§2. Five-dimensional quantum electromagnetism 

We consider the flat 5D space-time (X M ) = (x^,y) with the periodicity in the 
extra space y, 

ds 2 = rj flu dx fJ 'dx L/ + dy 2 , — oo < y < oo , y —* y + 21, 
(ry^) = diag(-l, 1, 1, 1) , (X M ) = (x»,x 5 = y) ee (x, y) , 

M,N = 0,1, 2, 3, 5; /x, i/ = 0, 1, 2, 3. (2-1) 

The 5D electromagnetism is described by the 5D U(l) gauge field Am, 

Sem = f d 4 xdyV^G{-\F MN F MN } = J d 4 xdyC EM , G = det G MN , 

F MN = d M A N -d N A M , (X M ) = (x»,y) , 
ds 2 = G MN dX M dX N , (Gmn) = diag(-l, 1, 1, 1, 102-2) 

It has U(l) gauge symmetry, 

A M ^ A M + d M A , (2-3) 

where A(X) is the 5D gauge parameter. 

We respect Z<i symmetry in the extra space, 

y^-y . (2-4) 

The Zi parity assignment of Am(x^, y) is fixed by the 5D gauge transformation (|2-3p . 
There are two cases corresponding to the choice of the Z2 parity of A(x fJ, ,y), 

Case 1. Even-parity case A(x tl ,y) = +A(x fl , —y) 
\ : P = + , A 5 : P = - , 
Case 2. Odd-parity case A(x^,y) = —A(x^, -y) 

A,j, : P = - , A 5 : P = + . (2-5) 

In the present paper, we consider Case 1. (Case 2 can be similarly treated.) 
We take the following gauge-fixing term to quantize the present system. 

C g = - l -{d M A M ) 2 = - l -{d^ + d y A") 2 , 

£em +£g = \ A ^ d2 + d y) A>M + \ A5 ( d2 + d y) A5 + total derivatives , (2-6) 
where d 2 = d^d^. Then the field equations are given by 

{d 2 + 8 2 )A^ = , (d 2 + d 2 )A 5 = . (2-7) 

*' The gauge independence of the physical quantities is an important check point of the proposal 
in the present paper. We relegate it to a future work. The gauge independence of Casimir energy 
of the 5D KK theory (Appelquist and Chodos's resu lt 3 >) was confirmed in Ref. ' 
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We consider the system in the periodic condition (|2-ip . Then we can write A M as 

oo _ 

A*{x,y) = <(x) + 2^a£(x)cos^y , P=+ 

n=l 1 

oo 

A 5 (x,y) =2^6 n (x)sin^y , P= - , (2-8) 

n=l 1 

where {an(x)} and {b n (x)} are the KK-expansion coefficients. From Eq. (|2-7p . they 
satisfy 



d 2 a% = (zero mode), {d 2 - (-y-) 2 }< = 0, {d 2 - {-j-) 2 }b n = 0, n ^ .(2-9) 

The on-shell condition, Eq. (|2-7p or Eq. (l2-9D . is for the analysis of the S-matrix. In 
this study, we do not use the condition. q| The total action can be written as 




*' We do not take into account the degree of freedom 5—2=3 among five components {^4 M } 
due to the local gauge symmetry. This is because we will compare the present results of the flat 
geometry with those of the warped geometry. In the latter treatment, we start with the massive 
vector theory that has no local gauge symmetry. 9 ^ 



S 5dv = J d 4 xdz^/—G(— —Fm n F M N - - m 2 A M A M ) , . (2-10) 

The 5D mass parameter m is regarded as an IR-regularization parameter. The 5D gauge theory 
is the limit m = 0. For the general m, Casimir energy is some integral of the (modified) Bessel 
functions with the number v = yT+ m 2 /oj 2 where uj is the warp parameter. The simplest case 
for analysis is not m = but m = iu). UV behaviour, which is the key point of the present 
paper, is considered independent of the IR regularization parameter m. Hence, the massive case is 
(practically) important for the 5D theories. We also include the longitudinal component to respect 
the manifest 5D Lorentz invariance. For the later use of the comparison with the warped case, we 
consider, instead of the 5D EM, the system of four 5D massless scalars with the even parity and 
one with the odd parity. 

**' Casimir energy is defined to be the free part (independent of the coupling) of vacuum energy 
that depends on the boundary. The quantity is defined to be the energy per unit space-volume of 
the "brane". In the present model of 3-brane (3+1 dim real world), space- volume has the dimension 
of I/ 3 . Hence the dimension of Ecas is L~ 4 . 
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where p 2 = p^p^ and m n = ^j-. This expression is the standard one. The above 
KK-summation and p^ integral are divergent, hence we must regularize them. The 
standard way taken by Appelquist and Chodos 3 ^ goes as follows. It is sufficient to 
consider the even-parity quantity. 



m r . 



(2-13) 



This is the (unregularized) Casimir energy for one scalar mode with Z2-parity even. 
The first step is to introduce a reference point ^o- 



V(l) - V{1 ) 



d 4 p 



(p3 + W) 



2 J (2vr)4A^ (p 2 + (^) 2 ) • 



(2-14) 



lo 



This procedure makes us drop the /-independent quantity. Using the well-known 
formula 



E /» 



dz f(z) + 



+co+le d J(z) + f(-z) 



-2iriz 



(2-15) 



the KK-sum in Eq. (|2T4p is replaced by the z-integral: 



v(i) - v(io) 



1 f d 4 p 



2 J (2tt) 4 



dz In — ' + 

P 2 + (f) 2 



+oo+ie 2 In 



p 2 + (?) 2 



(2-16) 



We consider the spacelike 4D momentum p„: p 2 = p^p^ > 0. Using the second 
formula, 



f°° z 2 + a 2 f a 

/ dz H(z) In 2 = 2?r / dx H{ix) , (2-17) 



z 2 + 6 2 

the first part of f)2- 16|) is evaluated as 
1 /• d 4 p 



2 J (2tt) 4 



°° dz l n ^±lxl!) 



(l - h) 



d 4 p 

(27T) 4 



(2-18) 



This integral is quintically divergent, but turns out to be cancelled out, as shown 
below. As for the second part, the integrand of the p-integral, using Eq. (l2-17p again, 
is evaluated as 



dz 



(p 2 +(^f) 2 ) 



-2-Ktz 



1 



2vr 



2nx 



-dx 
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We see that the first term in the above final expression, after the p-integral, cancel 
the quintically divergent one f)2- 18j) . Hence, we finally obtain 

V(l)-V(l )= In 



8vr 2 Z 4 J 
Using the third formula 



(2vr) 4 e l oy /^ _ e -ky/F 

/ dkk 3 {k + ln(l - e~ 2k ) ~Y k ~ ln (! ~ e ~~^)}> ( 1 Vp* = A;).(2-20) 
Jo ' 



f 

Jo 



we obtain 



1 3 

dk k 3 ln(l-e- 2k ) = -|C(5) , (2-21) 



^^Ff^^H^H^'-i-^} • (2 - 22) 

The first term is quintically divergent. We take, as the (dimensionless) UV cutoff of 
the /c-integral, M then 

m -V m = l^-^-^-\^) . ,2-23) 

Hence, we obtain Casimir energy and Casimir force for the simple system Eq. (|2-13p 



as 



5 4 / 4 ' ° asv 7 9/ V 5 I 5 J 8vr 2 ' 

C(5) = 1.03693 •• ■ .(2-24) 

The first term of V(l) is quintically divergent but is simply proportional to I. This 
quantity comes from the UV divergences of 5D quantum fluctuation. In Casimir force 
^Cas = ~ TJT> ^ na ^ P ar t does not depend on I. If we can find a right means of avoiding 
the UV divergences (which will be proposed later) we may drop the (divergent) 
constant contribution, and we obtain Casimir force for the 5D electromagnetism as 

F Oas = 5 x (-3^) x JL . (2 .25) 

The minus sign denotes the attractive force. 

§3. Heat-kernel approach and position/momentum propagator 

We reformulate the previous section using a heat kernel in order to treat the 
problem without KK expansion. Note that the heat-kernel method is a complete 
quantization procedure for the free theory (quadratic theory). 21 - 1 Instead of the 
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5D gauge fields A M (X), we introduce partially Fourier-transformed ones A¥ (y) 
(A$(y),B p (y)). 

d 4 p 



A»(x,y) = J ^e^A^y) : P=+ 



A 5 (x,y) = J ^e^B p (y) : P= - . (3-1) 

(We do not Fourier-transform the extra space (y) part.) Then the total action is 
given by 



ay 



(27T) 



S = j d 4 xdy(C E M + C g ) 
U^yX-p 2 + d y 2 )A^y) + \B p {y){-V 2 + d y 2 )B p (y) 



.(3-2) 



Here we restrict the y-integral region to [—1,1] because it has sufficient information 
and can be transformed (by the Fourier expansion) to the periodic form defined in 
(— oo,+oo). The on-shell condition is given by 

{-p 2 + d y 2 )A^(y) = , -l<y<l , P= + 
( -p 2 + d y 2 )B p {y) = , -l<y<l , P=- . (3-3) 

This condition, which is not used in the following, is necessary when we consider the 
S-matrix. Casimir energy E(j as is given by 

e -i 4 E Cas = J VA^VBpexpiiS} 

Using the following formula 21 ) 

dt = hxM , det M > , M : a matrix , (3-5) 



L 



oo e -t _ e -tM 



/0 t 
we can formally write 

- ln(p 2 - d y 2 ) = / -e- t( - p2 ~ d y ] dt + divergent constant , (3-6) 
Jo t 

where the divergent constant should not depend on p or y. We understand that M 
in Eq. ()3-5p is the matrix M^y labeled by the continuous parameters y and y' and 
that (p 2 — d y 2 ) in Eq. (|3-6p is the differential operator acting on \y >, a quantum 
state labeled by the position y. E3 The heat kernels H p and E p are in an abstract 

*' < y\ and \y > are introduced by Dirac 22 ' and are called the bra and ket vectors respectively. It 
is defined by the orthonormal eigenfunctions of the kinetic differential operator of the system. When 
we take the orthogonality relation < y\y' >= 5(y — y 1 ), their physical dimensions are [ — y^]=L _ 2 
and [jy — }—L~i . 
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way defined by 



H p (y,y';t)=<y\e-^- d >\y'> 
E p (y,y';t)=<y\e-^- d >\y'> 



(3-7) 



Hence, we obtain the final expression of Ecas'- 



(const) x 



cxp 



d^ P r°° dt f 4 



TV £p(y> y; t) + -Tr flp(j/, y'; t) 



(3-8) 



where Tr represents the integral over all y = y' values. E3 

TrE p (y,y';t) = J dy E p (y,y;t), Tr H p (y, y';t) = J dy H p (y,y;t). (3-11) 



The precise definitions of H p and E p , with the initial condition (|3-14p shown 
below, are given by the heat equations, 



{§- t +p 2 - d y} H p(y>y'> t ) = ° > p = - > 

{^+P 2 -^W^';*) = , P = + . 



(3-12) 



The solutions are, in terms of the KK-eigen- functions, given by 



neZ 



t.J .2\, 1 

2 



} , 



neZ 



(3-13) 



where we use the dimensionalities of H p and E p read from (|3-8p ; [£^]=[//p]=L . 



*' For the 5D free scalar with Z2-parity even, Casimir energy is given by 



e - cas _ ( cons t) x exp 
and, for that with Z2-parity odd, 

e - ; E cas _ ( C onst) x exp 



J (2tt) 4 J t 2 pVy ' y ; 



r d^ P r-dti , 



t) 



(3-9) 



(3-10) 
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The above heat kernels satisfy the following b.c, 

lim Q H p (y,y';t) = i £ I{ e -*»W) - e -^(y+y') } 

neZ 

= \{Kv - y) -<% + y')} , 

t lim o ^(y, y ';t) = i £ I{ e ~*-W) + e -<M*V) } 

= \{hy-y) + 5(y + y')} , (3-14) 

where we have introduced <5(y) = A ^ ngZ e~ lknV . With this b.c, the heat equation 
(|3-12|) can rigorously define H p and E p . We here introduce the position/momentum 
propagators G^ as follows: 

G;(y )? /0 = / dt H p (y,y';t) = - £ - -{e"*^ > - e -*My4V) }> 

/*oo i 11 

G+(y,y') = / dt E p (y,y';t) = - £ _-j {e^^ ) + e - ife »^')}.(3-15) 
They satisfy the following differential equations of propagators: 

neZ 

^Wy-y')T <% + </)} , (3-16) 
Therefore, Casimir energy Ec as is, from Eqs, (|3-8p and (|3- 13|) . given by 

+ 2 I^ e -(^+P 2 )*I{l + e -2^nS/}l _ (3 . 17) 

neZ J 

This expression leads to the same treatment as that in the previous section. 

Here, we introduce the generalized P/M propagators, J a (P=— ) and J a (P=+) as 



ia{p 2 ;y,y') = —H p (y,y';t) 

Jo 1 



**' If we ignore the dimensionality, 

j J2{e~ {k " +p2)t /(ki + p 2 ) s }^{e- tkn(y - y ' ) T e- lkrl(v+y,) }, s : real number 

ngZ 

are the more general solutions of Eq. (I3-12|) . 
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00 dt 1 
n ¥21 



a 07 9 l 



*>} 



neZ 



J a (p 2 ;y,y') 







—E p (y,y;t) 



00 1 



nez 



+ e -ifc„(2/+2/')} ^ p =+ ; (3 . 18 ) 



where a is an arbitrary real number. Then we have the following relations: 

h{p 2 \y,y') = G p (y,y) , J Q {p 2 ;y,y') = G p v (y,y') , 

dl a (p 2 ',y,y') r / 2 /\ l°° ji2 t n 2 /\ r ^2 /\ 
^2 = -la-liP ;y,y ), / dk I a (k ;y,y ) = I a+1 (p ;y,y ), 

dJ a (p 2 ;y,y') 



dp 2 



v 

00 



(p 2 - 9^)1 pip 2 ;y,y') = -Pi/3+i(p 2 ;y,y') , 

(P 2 ~ d 2 )Mp 2 ;y,y') = -PJ/3 +1 (p 2 ;y,y') ,(3^0 . (3-19) 



Finally, we obtain the following useful expression of Casimir energy in terms of P /M 
propagators, 



Ecas{l) 



d A p fl 



(2tt) 4 [2 



Tr h(p 2 ;y,y) + ^Tr Ji(p 2 ;y,y')| 
= f -0^ £ ^TrI (k 2 ;y,y')+2TrJ (k 2 ;y,y')}dk 2 

Here, we list the dimensions of the various quantities appeared above. 



(3-20) 



Dim 


L l-2a 


L- 4 


L -3/2 


L- 1 


L-1/2 


L 


L 2 


L 5/2 




la, J a 


Ecas 


A* 1 , A 5 


fi,A,p 
Hp, Ep 

s(y-y') 


\y >,<y\ 


u p 


t 


Ap , Bp 



(A, jj,, and e are the regularization parameters defined below.) 

The P/M propagators Gf, which are expressed in Eq. (|3-15p in the form of a 
summation over all modes, can be expressed in a closed form. (See, for example, 
Eq.(36) of Ref.10).) 



cosh k(\y + y'\ — I) T cosh k(\y — y'\ — I) 



4k sinh kl 

-I <y <l ,-l <y <l , k = ^Jk^kv , k^ > O(spacelike) , (3-21) 
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where the plural sign means that one corresponds to the other in the same position. 
By using the above results, Casimir energy is explicitly written as 

E C as(l) = l^iy y (F-(p, y )+AF + (p,y)) , 

cosh£;(2y — I) — cosh kl 



F~(p,y)^ I dk 2 G~ k (y,y) = I dk 



p 



2sinh(fc/) 



tt+,~ \ r° Jl2 ^+, \ f 00 ~- cosh k{2y -I) -cosh kl , n nn . 

This is the closed expression, not the expanded one. (In relation to the degree of 
freedom of the system, the integrands of F T are, in a later description, graphically 
shown in Figs JS] and [6l) 

§4. UV and IR regularization parameters and evaluation of Casimir 

energy 

The integral region of Eq. (|3-22p is displayed in Fig. 1. In the figure, we introduce 
the UV and IR regularization cutoffs ji < p < A and e<y <l. In order to suppress 
the number of artificial parameters as much as possible, we take the relations lj : 

e=i , *=\ ■ (4-D 

This is the same situation as that in the lattice gauge theory ( unit lattice size = 
(1/yl) 4 x e = (1/M) 5 ; total lattice size =(l//x) 4 x I = I 5 ). 
Let us evaluate the (A, l)-regularized value of Eq. (|3-22p . 

2^2 j-A r-l 

E Cas (A,l) = —— dp dyp 3 F(p,y) , 
( 27r ) Jl/l hi A 

% /~ n n f ,r — 3 cosh k(2y — I) — 5 cosh kl . 

F(p,y) = F-(p,y)+AF + (p,y)= / dk \ ' . (4-2) 

Jn 2smh (kl) 



The integral region of (p, y) is the rectangle shown in Fig. 1. We now use the formula 
of the indefinite integrals: 



/cosh, x 
— dx = In sinh(x) 
sinhx 



The expression of Eq. (|3-21|) is considered in — I < y < I, and the periodicity (y — > y + 21) 
seems lost. If, however, we Fourier- expand Eq. (|3-21l) in the interval, the same expression as Eq. (|3-15|) 
is obtained. Note that the present treatment is crucially different from the deconstruction approach 
in that all Kaluza-Klein modes are taken into account. The necessity of all the KK modes was 
stressed in the 5(0)-problem of the bulk-boundary theory. 23 ^' 24 ' 

*' As for the numerical values of A and I, they depend on the chosen energy unit ( e.g. eV 
and J) and the physical model concerned. For example, if we consider the grand unified theories, 
A = 10 19 GeV (Planck energy) and r 1 = 10 3 GeV (Tev physics) are strong candidates. Another 
case is the cosmological model. Then we take A = 10 19 GeV and r 1 = 1(T 41 GeV ([Cosmological 
Size] -1 ). A and I are huge numerical numbers when applied in the real world. 
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cosh(ax) 
sinhx 



(.1- a)x 



1 



-2Fl 



a 3 

1 J- 1 



-;e 



-2x 



-(l+a)x 
1 + a 



"2-^1 



1 + a ., 3 + a 



2 

x > 



1 



-;e 



dx 



-2x 



1 > a > , (4-3) 



where 2-^1 (a, /3; 7; 2) is Gauss's hypergeometric function. The formula below goes to 
the above one in the limit a — > 1 — (ignoring the divergent constant l/(a — 1)). 
The integrand of Ec a s(A,l) (j4-2j) . p s F(p,y), can be exactly evaluated as 



p 3 



fF(p,y) = -^- 



3 e - ^ 1-0 ^)) 1 
2 l-a(y) ; 



1 



o ' 5 



-2x 



3e 



-{l+a(y))x 



2 l + a{y) 



1F1 



1 + a(y) 3 + o(y) _ 2x 



+ — In sinh(x) 



x=Al 



I 



A' 1 <y <l 



x=pl 

r 1 <p<a 



(4-4) 



Note that Eq. ()4-2p . with Eq. ()4-4p . is the rigorous expression of (A, ^-regularized 
Casimir energy. We show the behaviour of p 3 F(p, y) by taking the unit I = 1 in 
Figs. 2-4. Three graphs correspond to A = 10, 100, and 1000. All three graphs have 
a common shape. The behaviour along p-axis does not markedly depend on y. A 
valley 'runs' parallel to the y-axis with the bottom line at a fixed ratio oip/A ~ 0.75. 
Rl The depth of the valley is proportional to A 4 . Because Ecas is the (p, y) 'flat- 



*' The valley-bottom line ('path') p = p(y) w Q.75A corresponds to the solution of the minimal 
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Fig. 2. Behaviour of p 3 F(p,y) (g3J. I = 1, A = 10, 0.1 < y < 1, 1 < p < 10 . 




Fig. 3. Behaviour of p 3 F(p, II = 1, yl = 100, 0.01 < y < 1, 1 < p < 100 . 

plane' integral of p 3 F(p, y) , the volume inside the valley is the quantity Ecas ■ It 
is easy to see that Ec a s is proportional to A 5 . This is consistent with Eq. ()2-24j) . 
Importantly, Eq. (|4-2p shows the scaling behaviour for large values of A and I. From 
a close numerical analysis of the (p, y)-integral Eq. (|4-2p . we confirm (see Appendix 



'action' principle: SS = 0, S= (1/87T 2 ) J dp J dyp 3 F(p,y). This will be referred to in $7] 
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Fig. 4. Behaviour of p 3 F(p,y) |g3}. 1 = 1, A = 1000, 0.001 < y < 1, 1<P< 1000 . 



C ) that 

27T 2 

EcasiAl) = [-0.1249M 5 - (1.41,0.706,0.353) x 10~ 5 M 5 ln(Zyl)] .(4-5) 

(Note: 0.125 = (l/8)x (1/5) x (4[ewen] + l[odd]).) This is the essentially the same 
result as that in ^21 l*)| 

Finally we notice, from Figs. [2HH that the approximate form of F(p, y) for large 
values of A and I is given by 

F{p,y)^-LA(l-^\ , f = 5 , (4-6) 

which does not depend on y or I. / is the degree of freedom. From this result, 
we know, using the integral expression (|3 22j) . that the following approximate form 
of the integrand of is valid in a wide range (k, y). 

, _ / ~ cosh k(2y — I) — cosh kl 1 
lntgrdF-(k,y;l) = \\ « -- , 



*' Note that the result of Eq. (|2-24p in Sec[2]is obtained from the KK-expanded form, whereas 
that of Eq. (|4-5|) in this section from the closed expression. The coincidence (A 5 proportionality) 
strongly shows the correctness of both evaluations. At the same time, the numerical result for the 
first term coefficient shows that the number of significant figures is 4. Altough the second term 
contibution is small, its coefficient is not stable and its significant digits is 1 at most. The triplet 
data correspond to I = 10, 20 and 40. 

*' Using the approximate form Eq. (|4-6I) . Ecas is estimated as Ec as (A,l) x 87r 2 w 
fi/i d P Ii/a d yf(- 5 / 2 )( A - P) = ~(1/8)M 5 (1 + O(^)). This is consistent with Eq.gSJ. 
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Intgr<±F + (fc, y; I) = 



cosh k(2y — I) — cosh kl 
2sinh(fci) 
(k,y) G > 1 and fc(Z - 



y) » 1} 



(4-7) 



which can be analytically proved and is confirmed by graphically showing the above 
functions (see Figsj5] and [6j The table shape of the graphs implies the "Rayleigh- 
Jeans" dominance becaus e C asimir energy density is proportional to the cubic power 
of p in the region p -C A. \**' 




100 



Fig. 5. Behaviour of the integrand of F~, IntgrdF ~ (k, y; I) (|4T|) . I = 1, A = 100, < y < I = 1, 
1 < k < A = 100 . The fiat plane locates at a height of —0.5. 



§5. UV and IR regularization surfaces, principle of minimal area and 

renormalization flow 

We have confirmed that heat-kernel formulation is equivalent to the familiar KK- 
expansion approach. The advantage of the new approach is that the KK expansion 
is replaced by the integral over the extradimensional coordinate y and all expressions 
are written in the closed (not expanded) form. The yl 5 -divergence shows notorious 
problem of higher dimensional theories. In spite of all previous studies, we have 
not succeeded in defining higher-dimensional theories. (In this free theory case, the 
divergence is simple and the finite Casimir energy (force) can be read as shown in 

**' The well-known radiation spectral formula (/3 is the inverse temperature): {E) Vl g = 
hv/2 + hvl{e^ hv — 1) consists of two parts. The first one is the zero-point energy and the 
second is Planck's ditribution part. The low-frequency region of Planck's dstribution formula: 
hv/(er hv — 1) w 1/(3 (independent of v), hv <<iC 1/(3, is called Rayleigh- Jeans region. Note that the 
extra coordinate y or I — y, in the present 5D model, plays the role of the inverse temperature (3. 
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Fig. 6. Behaviour of the integrand of F+, IntgrdF+(fc, y; I) l|¥7|l . I = 10, A = 10, < y < I = 10, 
1 < k < A = 10 . The fiat plane locates at a height of —0.5. 



Eq. (l2-25p . In general, however, the divergences cause problems. The famous example 
is the divergent cosmological constant in the gravity-involving theories. 3 ) ) We notice 
here that we can solve the divergence problem if we find a way to legitimately restrict 
the integral region in (p, y )-space. 

One proposal of this was presented by Randall and Schwartz. 9 ) They introduced 
the position- dependent cutoff, \i < p < 1/u , and u G [e,l] , for the 4D-momentum 
integral in the "brane" located at y = u. (See Fig. [TJ) The total integral region is 
the lower part of the hyperbolic curve p = 1/y. They succeeded in obtaining the 
finite /9-function in the 5D warped vector model. We have confirmed that the Ec as 
obtained using Eq. (|4-2p . when the Randall-Schwartz integral region (Fig. [T|) is taken, 
is proportional to A 4 . A close numerical analysis shows (see Appendix C) 

27T 2 /"V"? 27T 2 f 1 /" 

E c!s = 7^ dq dy q 3 F(q,y) = -— [ du dpp 3 F(p,u) 
Ji/i hi a ( 27r ) hi a Ji/i 

2tt 2 

= [-8-93814 xlO- 2 ^ 4 ] ,(5-1) 

which is independent of I. This shows that the divergence situation is indeed 
improved compared with the yl 5 -divergence of Eq. (|4-5p . 

Although Randall and Schwartz claim that holography (for the case of the 
warped geometry) is behind the procedure, the legitimateness of the restriction looks 

*' The result of Eq. (|5-1[) is consistent with the approximate form of F obtained in Eq, (14-6[) . 
(-5/2) dqfljl dyq 3 (A - q) = -(l/12)yl 4 (l + 0(l/(yl/) 3 )). 0.0893 » 0.0833 • •• = 1/12. 
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less obvious. We have proposed an alternative approach and given a legitimate expla- 
nation within the 5D QFT. 10 )' 25 )' 26 ) Here we closely examine the new regularization. 



P t 



A -■ 



[X - 




Y 



Fig. 7. Space of (p,y) for the integration (present proposal). 



On the "3-brane" at y = e, we introduce the IR cutoff [i and the UV cutoff A 
(fi <C A) (see FigE]), 

M < A . (5-2) 

This is legitimate in the sense that we generally perform this procedure in 4D renor- 
malizable thoeries. (Here we are considering those 5D theories that are renormaliz- 
able in 3-branes. Examples are 5D electromagnetism (the present model), the 5D 
<£ 4 -theory, and the 5D Yang-Mills theory.) For the same reason, on the "3-brane" 
at y = I, we may have another set of IR and UV cutoffs, // and A' . We consider the 
case, 

// < A', A' < A, n ~ n' . (5-3) 

This case will enable us to introduce renormalization flow, (see our later discussion.) 
We claim here that, as for the "3-brane" located at each point y (e < y < lY 
the regularization parameters are determined by the minimal area principle. q| 
To explain this, we move to the 5D coordinate space (see Fig. [8]). The 

p-expression is replaced with the ^/x^rr^-expression by reciprocal relation, 



\Jx ll {y)x^{y) = r(y) 



p(y) 



(5-4) 



*' We do not quantize the (bulk) geometry, but treat it as the background. The (bulk) geometry 
fixes the behavior of the regularization cutoff parameters in the field quantization. The geometry 
influences the "boundary" of the field-quantization procedure in this way. 



18 



S. Ichinose 




The UV and IR cutoffs change along the y-axis and the trajectories form surfaces 
in the 5D bulk space We require that the two surfaces do not cross for the 

purpose of the renormalization group interpretation (discussed later). We call them 
UV and IR regularization (or boundary) surfaces {Bjjy and Bjr), 



1 



B UV : vVp + (x 2 ) 2 + (x 3 ) 2 + ( X 4 )2 = ruv{y) ; e = _ <y</ , 

B/H : ^(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 + (x 4 ) 2 = r IR {y) , e = ^<y<Z ,(5-5) 

The cross sections of the regularization surfaces at y are the spheres S 3 with the 
radii ruy(y) and rm(y). Here, we consider the Euclidean space for simplicity. The 
UV-surface is stereographically shown in Fig. [S] and reminds us of the closed string 
propagation. Note that the boundary surface B[/y (and B/#) is the four-dimensional 
manifold. 

The 5D volume region bounded by Bjjy and Bir is the integral region of Casimir 
energy Ecas- The forms of rjjv{y) an d rj^(y) can be determined by the minimal 
area principle, 

^(Surface Area) = , 3 % — = , < y < I . (5-6) 

In Appendix A, we present the classification of all solutions (paths) and the general 
analytic solution. In Fig. [101 we show two result curves of Eq. (15-6p . taking the 
boundary conditions (r' = dr/dy): 

Fig. QUI Fine configuration goes to coarse configuration as y increases 
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IR curve (upper): r[0] = 12.0, r'[0] = —1.0 [simply decreasing type], , , 
UV curve (lower): r[1.0] = 10.0, r' [1.0] = 350.0 [simply-increasing type]. 

where the types of curves are specified on the basis of the classification of the minimal 
surface lines in App. Al. They show that the renormalization flow shown in Fig. [7] 
really occurs by the minimal area principle. In Fig. llll another set of minimal surface 
lines are given by taking another boundary conditions. 

Fig. [TTJ Coarse Conf. goes to Fine Conf. as y increases 
IR curve (upper): r[0] = 4.6, r'[0] = —1.0 [simply decreasing type] , , . 
UV curve (lower): r[0] = 4.5,r'[0] = —22.0 [simply decreasing type] 

They show the opposite-direction flow of renormalization compared with that in 
Fig. Q3U (See the next paragraph for the renormalization flow interpretation.) These 
two examples imply that the boundary conditions determine the property of the 
renormalization flow. E3 

The present regularization scheme also gives the renormalization group interpre- 



tation of the change in physical quantities along the extra axis (see Fig. [8j). t—1 In 



*' The minimal area equation Eq. (|5-6I) is the second-derivative differential equation and has 
the general solution (see Appendix A). Hence, for the given two initial conditions (for example, 
r(y — e) and dr/dy\ y=e ), there exists a unique analytic solution. The presented graphs are those 
with these initial conditions. Another choice of the initial conditions, r(y = e) and r(y = I), 
is possible. Generally, the solution of the second-derivative differential equation is fixed by two 
boundary conditions. 

**' This part is in contrast to the AdS/CFT approach where the renormalization flow comes 
from the Einstein equation of 5D supergravity. 





1 



0.2 0.4 0.6 0.8 



Fig. 11. Numerical solution of (|5-6[) . Vertical axis: r; Horizontal axis: < y < I = 1. IR curve 
(upper): r[0] — 4.6, r'[0] = —1.0 [simply decreasing type]; UV curve (lower): r[0] = 4.5, r'[0] = 
—22.0 [simply decreasing type]. 
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the "3-brane" located at y, the UV cutoff is rjjy{y) and the regularization surface is 
the sphere S 3 with the radius Wviv)- The IR cutoff is rjn(y) and the regularization 
surface is another sphere S 3 with a radius rm(y). We can regard the regularization 
integral region as the sphere lattice of the following properties: 

Unit lattice (cell): the sphere S 3 with radius ruv(y) an d its inside, 
Total lattice: the sphere S" 3 with radius rjn(y) and its inside. 

It is composed of many cells above, 

(TjjAy) \ 4 
tt • (5'9) 
ruvWj 

The total number of cells changes from (^) 4 at y = e to (^r) 4 at y = I. Along the 
y-axis, the number increases as 

rm(y) i - N ( y ) . (540) 



juviy) 

For the "scale" change y — > y + Ay, N changes as 



When the system has some coupling g(y), the renormalization group (3{g) function 
(along the extra axis) is expressed as 

A(lng) 1 Ag _ 1 1 1 dg 

P ~ A{\nN) ~ A(lnN)T~^9 H nRhj) ) gdy' ' (542) 

ay K r uv (y)> 

where g(y) is the renormahzed coupling at y. 

We have confirmed that the minimal area principle determines the flow of the 
regularization surfaces. Among the numerical results, some curves are similar to the 
type proposed by Randall- Schwartz. 

§6. Weight function and Casimir energy evaluation 

In Eq. (|3-22j) . Casimir energy is written by the integral in the (p, y) space in the 
range: < y < I and < p < oo. In § [SJ we see that the integral region should be 
properly restricted because the cutoff region in the 4D world changes along the extra 
axis obeying the bulk (flat) geometry (minimal area principle). In this section, we 
consider an alternative version. 

Instead of restricting the integral region, we introduce the weight function W(p, y) 
in the (p, y)-space to suppress the UV and IR divergences of Casimir energy. 

/»4 pi 
—^4 J q dyW(p,y)F(p,y), F(p, y) = F~ (p, y) + AF + (p, y) , 

*' Here we consider interacting theories, such as the 5D Yang-Mills theory and 5D <P 4 theory, 
where the coupling g(y) is the renormahzed one in the '3-brane' at y. 
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Examples of W(p, y) : W(p, y) 
j_ 



( ^ e -(i/2)(¥-(V% 2 /i 2 = Wl (p,y), JVi = ^ 

820 



^e-™ = W2(p,y), iV 2 



2(L4) j 



0.3222 



iV 6 

i e -(l/2)lV 



Sir- 



T7 , f W 7 (p, y ), AT 7= ™ 



3800 



fe-Q'/iHP+i/v*) = W 8 {p,y), N 8 = -^ 
_i_ e -(/W+i/, 2 ) 2 = ^ 88fey) , at 88 = a— 

03921 



elliptic sup., 
kinetic-energy sup., 
hyperbolic sup. 1, 
hyperbolic sup. 2, 

linear sup., 
parabolic sup. 1, 
parabolic sup. 2, 
higher-derivative sup. 1 

reciprocal sup. 1, 
higher-derivative sup. 2 

reciprocal sup. 2, 
higher-der reciprocal sup 
reciprocal sup. 3. 



where are defined in Eq. (|3-22p and iV, are the normalization constants. — 



Above, we list some examples expected for the weight function W(p,y). W% and W3 
are considered to correspond to the regularization taken by Randall-Schwartz. How 
to specify the form of W is the subject of the next section. We show the shape of the 



*' In the flat geometry, the periodicity parameter I is a unique scale parameter. We make all 
exponents in Eq. (|6-1[) dimensionless using I. 

**' The normalization constants N are defined by 

f rl^yi 1 1 pi A pi i 

/ 77^1 / d y Wi &^ = ^7? / dx / dw * a Wi(y.H = 73 - lA » 1 ( 6 '- 2 ) 



where x = Ip and w = y/l. They are explicitly given by 



Stt 2 ^ 



dw e~ w2/2 = 1.557 , 8tt 2 7V 



87r^JV 2 



00 pi 



d.i: 



(iA)-i 



dwx 3 e~ xw = 2(M) 3 , 87r 2 iY 3 



87T 2 iV4 



87T 2 iV 8 ' 



Ve Jo 



dw e~ 1/2w2 = 0.38 00 , 87r 2 iY 4 - 



871-^56 



8nN 88 



87r 2 7V lb 


3 f°° , 3 
= — = = / dx x e 






/*oo 




=y i *, 


f 1 j 3 -i 2 iu 2 

/ dwx e 




•/(M)-i 




/■oo p 1 


87T 2 iV 5 = 


/ da; / dw x 3 e~ 








f'°° 4 


87r 2 7V 7 = 


/ dx x 3 e~ x ^ 2 — 








1 


f dx I 


j 3 -x 2 lx 2 + l/w' 

dw x e 



3 -x 2 /2 



1.820 , 



(2/3)(M) 3 , 



2^ 



0.3033 



dx / dw a;3 e -(-/»)(-+i/™)/2 = al346 



^ dx f 1 dw s 3 e -(l/2)(x 2 +l/^) 2 = 0i005006 



00 /•! 



87r 2 7Vg = / dx dw a; 3 e -( 1 /2)(^+i/») = 0.03921 .(6-3) 
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energy integrand p s W(p, y)F(p, y) in Figs. [T2TU51 We notice that the valley-bottom 
line p ~ 0.75/1, which appeared in the unweighted case (Figs .2-4), is replaced by new 
lines: p « const (Fig. [12] ,Wi), py ~ const (Fig. [HIW3), p ~ const x y (Fig. EBW4), 
p « const x y/y (Fig. dSJWg). They are all located away from the original /1-effected 
line: p ~ 0.75/1. 




Fig. 12. Behaviour of p 3 Wj (p, y)F(p, y)(elliptic suppression), yl = 10, I = 1, 1/yl < y < 
0.99999/, l/l <p< A. 



We can check the divergence (scaling) behaviour of E^ as by numerically evalu- 
ating the (p, y)-integral (|6-ip for the rectangle region of Fig. 1. □! 



See Appendix C for the explanation of the numerical derivation. 
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Fig. 13. Behaviour of p 3 Wz{p,y)F{p,y) (hyperbolic suppression 2). A — 10, / = 1, 1/A < y < 
0.99999/, l/l <p< A. 



-(2.500, 2.501, 2.501) f| + (-0.142,1.09,1.13) x 10^ 4 /11 ^ M) for W u 

-(2.502, 2.502, 2.502)^ + (2.40, 2.44, 1.84) x i - 4 il^l f or W lb , 
-(6.0392,6.0394,6.03945) x 10™ 2 ^ 

-(24.7,2.79,1.60) x 10~ 8 ^^) for W 2 , 
-(10.650,9.21917,9.21915) x 10~ 2 ^ 

+(153.3,1.9629624,1.9629620) x W~ 5 Aln ^ A) for W 3 , 

-(2.55, 2.53, 2.52)f| + (4.56, 2.63, 1.48) x io-s^M^) for ^ 

-(2.55,2.54,2.52)^ + (6.10,3.63,2.10) x 10^ 3 ^||^ for W 5 , (6-4) 

-(2.532, 2.519, 2.511)^ + (3.19, 1.83, 1.03) x ltr 3 j^*) f or w 6 , 

-(2.51, 2.51, 2.50)$ + (8.51, 5.51, 3.36) x 10-* Ah $ A) for w ^ 

-(2.52, 2.51, 2.51)$ + (19.5, 11.6, 6.68) x lO" 4 ^^ fo r w 8 , 

-(2.55, 2.55, 2.55)$ + (5.47, 5.32, 5.01) x iq-3 ^M^) for WiJf 

-(2.54, 2.53, 2.52)$ + (4.30, 2.48, 1.40) x 10~ 3 ^(^) f or w 56 , 

-(2.61, 2.61, 2.60)$ + (1.20, 1.17, 1.10) x W ~2Ah^A) for ^ 

-(2.54, 2.52, 2.51)4 + (3.75, 2.16, 1.22) x 10~ 3 ^^ for W 9 . 



The above fitting is obtained by taking the data for the range: I = (10,20,40), 
A = 10 ~ 10 3 . The round-bracketed triplet data, corresponding to three values of I, 
should be the same if the scaling region is properly examined. The small fluctuation 
in the last digit number tells us the significant figures. The leading terms (linearly 
divergent terms) are firmly obtained, and the number of the significant figures (NSF) 
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Fig. 14. Behaviour of p 3 W4(p, y)F(p, y) (linear suppression). A = 10, I = 0.5, 1/A < y < 
0.99999/, l/l <p< A. 



> 3. The log terms are obtained very poorly, and the NSF is 1 at best. The 
divergent nor maliz ation constants of W2 and W3, namely, N2 and N3, are consistent 



with Eq. (|5-ip . 1**1 After normalizing the factor Al, only the log-divergence remains. 



E^ as /Al = -^(l-4cHlA)) , (6-5) 

where a can be read from Eq, (|6-4p . (For the 5D EM (|2-24j) . a corresponds to the value 
5 x 3£(5)/(4 x 8ir 2 ) ~ 3.86/87T 2 .) The numerical results say that a does not so much 
depend on the choice of W except W2 and W3. As for c, we expect that it reaches a 
fixed value as I increases further. Although the present approach leaves the weight 
function W(p, y) unspecified, and the numerical results involve some ambiguity, we 
can say that a is approximately +2.5 (a positive number on the order of 0(1)), and c 
is a positive number on the order of O(10 -4 ) to O(10 -3 ). (The hyperbolic cases, W2 
and W3, are exceptions. This is due to the divergent normalization factors N2 and 
./V3 shown in Eq. (|6Tp ). At present, we cannot discriminate which weight is the right 
one. Here, we list the characteristic features (advantageous (Yes) or disadvantageous 
(No), independent (I) or dependent (D), and singular (S) or regular (R)) for each 
weight from the following points. 

point 1 The behavior of W for the limit I — > or / — > 00. This property is related 
to the continuity to the ordinary field quantization. 

point 2 The path (bottom line of the valley) is independent (I) of the scale I or 

*' No data instability appears in a warped case. 26 - 1-28 - 1 
**' In particular, Wi case is more similar to that in Eq. (|5-1[) in that the log term vanishes in the 
present numerical precision. 
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Fig. 15. Behaviour of p 3 We(p, y)F (p, y) (parabolic suppression 2). A — 10, I = 0.5, 1.001/vl < 
y < 0.99999Z, l/l < p < A. The contour of this graph is given later in Fig. [201 



dependent (D) on it. 
point 3 Regular (R) or singular (S) at y = 0. 
point 4 Symmetric for Ip <-> y/l. 

point 5 Symmetric for p *-* 1/y, (Reciprocal symmetry), 
point 6 The Casimir energy is finite, 
point 7 Value of a. 

point 8 Under Z2-parity: y <-> — y; W(p, y) is even (E), odd (O) or none (N). 
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So far as the legitimate reason of the introduction of W(p, y) is not clear, we 
should regard this procedure as a regularization for defining higher-dimensional theo- 
ries. We give a definition of W(p, y) and a legitimate explanation in the next section. 
It should be conducted, in principle, consistently with the bulk geometry and the 
gauge principle. 



§7. Definition of weight function and dominant path 

In the previous section, the weight function W(p, y) is introduced as some trial 
functions for suppressing the UV and IR divergences. In this section, we define (or 
specify) the weight function W(p, y) properly and give a legitimate reason for the 
introduction of W. 

First, the requirement for controlling the yl 5 -divergence in § [J] led us to intro- 
duce some "damping" function W(p, y), as in Eq. (|6-l|) . Casimir energy is obtained by 
integrating out W(p, y)F(p, y) over the entire(5D) space region. Among all configu- 
rations involving the integral, there exists a dominant configuration that contributes 
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to the integral most dominantly. The present claim is that the dominant configura- 
tion should be fixed by the 5D geometry, that is, the 5D flat space(-time) with the 
periodic boundary condition and ^-symmetry. In the integral Eq. (|6-ip . the dom- 
inant path p = pw(y) is characterized by the differential equation obtained by the 
variation method: p — > p + 5p, y — > y + 5y in the (p, y)-integral expression Eq. (|6-ip . 
Here, we require, on the basis of the present claim, that the dominant path pw(y) 
coincides with the curve r = r g (y) (or its momentum counterpart p g (y), (|B-4j) ) which 
is determined by the minimal area condition E g. (15-61) . Note that the minimal area 
curve r = r g (y) is defined by the 5D geometry, □j 

To explain the previous paragraph using only the coordinate (r = y/x a x a , y) , we 
move to the coordinate expression by partial Fourier transformation. Casimir energy 
in Eq. (|6-ip is re-expressed as 

F(r(x),y) = J -0^e^F(p,y) , W(r(x),y) = J -0^e^W(p,y) , 

r(x) = ^(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 + (x 4 ) 2 , 
E^asil) = I d 4 x f dy W(r(x),y)F(r(x),y) = 



o 



rl poo 

2ir 2 dy dr exp{3 In r + In W(r,y) + In F(r,y)} .(7-1) 
Jo Jo 



The unweighted case is W(r(x),y) = S 4 (x), W(p, y) = 1. The dominant contribution 
(path) r = rw(y) to E^ as is given by the minimal 'action' principle, that is, by 
applying the steepest-descend method to Eq. (|7-ip . 

_J_dW_ _ 1 dF 

dr = wdy^ pay ^ y l{ w,d r W,d y W ;r,y) . (7-2) 

r y/ dr ' p dr 

(The valley-bottom lines that appeared in Figs. fT2Tll5l are regarded as the dominant 
paths.) Using the above result, we can obtain d 2 r/dy 2 . 

1 dW l dF - 1 dW I OF 
W dy F oy W or F or 

d 2 r _ dyW y Wy(d r W y + 8 y W r ) W 2 (-3r~ 2 + d r Wr) 



dy 2 3 r -i + w r (3 r -i + yy r )2 ( 3r -i + yy r )3 

= V 2 (W,d r W,d y W,d 2 W,d y d r W,d 2 W-r,y) . (7-3) 



*' When W(p, y) = 1, the dominant path appears as p{y) ~ 0.75yl (ind. of y) in Figs. 2-4. This 
path, however, is "artificially" created by the UV cutoff. It is irrelevant to the 5D geometry. The 
differential equation of pw{y) is obtained in Eq. (lB-6[) . 

**' In § 5, we required the minimal area condition on the UV and IR regularization surfaces 
(boundary configuration), whereas, in this section, we require it on the dominant configuration in 
the (p, y) or (r, j/)-integral appeared in the expression of Eq. (|6-f[) . In other words, we have fixed 
the 'dominant' configuration (path), around which a small ("quantum") fluctuation may occur, by 
taking the minimal surface curve. 



Casimir Energy of 5D Electromagnetism 



29 



We require here that the path r = rjy(y) of (|7-2p . which is defined in a W^r, y)- 
dependent way, coincides with the minimal surface curve r g (y) (15-6p . which is defined 
independently of W(r,y). Hence, W(r,y) is defined by inserting Eqs. (|7-2|) and (|7-3p 
in Eq. (15311: 

V 2 (#, d r ^, <9 r 2 #, d y <9 r #, d 2 #; r, y) 

--{Vi(W' 1 WAW , ;r,y)} 2 -- = . (7-4) 
r r 

We call this equation the "W-defining equation". It defines the weight function 
W(r,y) in terms of the bulk metric (geometry) and model information F. (Note: 
F is given by Eq. (|4-2|) .) In Appendix B.2, we treat the W-defining equation in 
(p = ^p a p a , y) variables. Besides, how much the trial weight functions satisfy the 
above definition is numerically examined. 

The scaling of the renormalized coupling g{y) is given by 

g _ 1 1 Idg 

i9_\nr(y)gdy ' 

where g(y) is a renormalized coupling at y. In the above derivation, Eq. (|5T2p is 
used in the case §^r IR (y) = 0, and r uv (y) = r(y). 

§8. Discussion and conclusion 



Let us suppose that we have found the right weight function and that the diver- 
gences are successfully suppressed logarithmically. Casimir energy (I2-24D is replaced 
by 

8^c as = -^(l-4cln(M)) = -^ , (8-1) 

where c is some constant (see Eq. (|6-5|) ). This shows that the periodicity parameter 
(or the compactification size) I changes as the renormalization scale changes. The 
parameter I suffers from the renormalization effect. It shows the field's interaction 
with the boundaries. The above relation is exact because the present system is the 
free theory, and the heat-kernel (1-loop) approach can be regarded as a complete 
(nonperturbative) quantum treatment. 21 ) Note that, in familiar 4D renormalizable 
interacting theories, the 1-loop effect is proportional to (coupling) 2 . In the present 
case, however, c is a pure number. When c is regarded small, i.e., c <C 1, we can 
approximate I' asE3 

V « 1(1 + cln(lA)) , (8-2) 
The scaling behaviour of / is given by 

A (/3-function) = — ^— In ~ = c . (8-3) 
o(mA) I 

*> From the results of Eq.([6^3|, c ~ O(10" 3 ). 
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When c > 0, the compactification size I grows (shrinks) as the cutoff scale A increases 
(decreases), whereas when c < 0, I shrinks (grows) as the cutoff scale A increases 
(decreases). The former case is expected. When c = 0, it means that the size I has 
no quantum effect. 

The present 5D geometry is flat. The interesting application is the warped case. 
The analysis is under way. Partial interesting results are obtained. 26 )~ 28 ) More or 
less, the arguments go similarly to that in the flat case except that the periodic 
and hyperbolic functions are replaced by the Bessel and modified Bessel ones. The 
essential difference is that one additional massive parameter, the 5D AdS curvature, 
appears in the expressions. Ecas is expressed by the massive parameter in addition 
to I and A. 

In the present standpoint, the space-time geometric field Gmn is regarded as a 
background one. It is not quantized. As for other bulk fields, we assume that they are 
renormalizable in the 3-brane. The role of the geometry appears when it is required 
that the dominant 'path', determined by the (EM) field quantization, coincides with 
the geometrically determined 'path' (minimal area principle). Practically W plays 
the role of suppressing the integral by weighting the original integrand. Although 
we have already stated the definition of W in § [7] and Appendix B, it is important 
to know the true meaning of W . In the next paragraph, we argue one possible 
interpretation. 

In order to most naturally accomplish the above requirement, we can go to a 
new step. That is, we propose to replace the 5D space integral in Eq. (|7- 1|) with 
the following path integral. Namely, we newly define Casimir energy in the higher- 
dimensional theory as follows: 

£ C as(l,A) = 



i dp f X\Vp a (y)F(p,y) exp 

Ji/A Jm)=m)=i/otz 



ll/A Jp(0)=p(l)=l/p a>y 



f iJlr + 1 d V 



2a' J p 3 V p 4 



[ dp f T\Vx a (y)F(-,y) exp L-L I vV 2 + 1 r 3 dy 

Jl/A 7r(0)=r(I)=P £,» r L 2 « io 



, (8-4) 



where the limit Al — ► oo is taken and the surface (string) tension parameter l/2a' is 
introduced (note: the dimension of a' is [length] 4 ). F(p,y) or F(l/r,y) is the energy 
density operator induced from the quantization of 5D EM fields in Eq. (|3-22p . The 
weight factor comes from the area suppression: 

exp(— Area/2a') = exp[— (l/2a') / ^detg a bd 4 x]. In the above expression, we have 
followed the path-integral formulation of the density matrix (see Feynman's text 29 -*). 
The above definition clearly shows that the 4D space coordinates x a or the 4D 
momentum coordinates p a are statistically qu ant ized with the Euclidean time y and 



the "area Hamiltonian" H = J \J det g a b d 4 x. fcj This reminds us of the space-time 



*' In the unified theories based on the string theory, I is called the moduli parameter and is 
given by the vacuum expectation value of the dilaton field. How to fix the parameter is the moduli 
stabilization problem. 

**' The possibility of the quantum feature of the 5D coordinate/momenta was pointed out in 
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uncertainty principle^ introduced in the development of the string theory. The 5D 
quantum field theory leads to some quantum statistical system of the 4D coordinates 
{x a (y)} with the inverse temperature parameter y. In this view, the treatment of §[6] 
and [7] is an effective action approach. We expect the direct evaluation of Eq. (|8-4p . 
numerically or analytically, leads to a similar result. 

In the convincing regularization procedure, there should be some limiting pro- 
cess. In the explanation of the proposed one, before the previous paragraph, we do 
not clearly state the limiting process. We take large, but finite, values of A and I 
and obtain some formulae through numerical calculation valid for general (large but 
finite) A and I. The limiting process should be given by 

lA — ► cxd , (8-5) 

as in the proposed Eq. (|8-4p . In fact this limit is taken in the evaluation of the nor- 
malization constants of the weight functions in Eq. (|6-ip . When the limiting process 
Eq, (|8-5p is taken into account, the physical quantity, that is, Casimir energy Ec as 
(Eq. (|8-ip ) has the log divergence ln(Zvl). This divergence is subtracted by the renor- 
malization procedure (|8-ip and (|8-2p . (Note that the system is a locally free theory. 
It does not have the mass and coupling parameters. The original divergences must 
be absorbed away by the renormalization of the wave function or boundary param- 
eter (s).) The anomalous scaling behavior Eq. (|8-3p shows the non trivial interaction 
between the (free) fields and the boundaries. 

In the present analysis, from the beginning, the extra space is treated differently 
from the 4D real space(-time). It is regarded as the part that provides the axis for 
a scale change. The change of renormalization group is determined by the minimal 
area condition. In order to extend the quantum field theory, without the divergence 
problem, the extra space should play a role in suppressing the singular behaviour. 
Although the geometry is treated as a background, the final outcome looks to demand 
some type of quantization among the 4D coordinates (or momenta) as explained in 
a previous paragraph. The necessity of the weight function W(p, y) in Eq. (|6-ip can 
be interpreted to mean that we must take the well-defined space (-time) measure 
d A p dy W(p,y), instead of d A p dy (Eq. (|3-22p or (|4-2p ). in summing over 5D space(- 
time). 

We have focused only on the vacuum energy. We must, of course, examine other 
physical quantities such as S-matrix amplitude in the interacting theories. 

The present proposal should be compared with the string theory. We do not 
directly treat the string propagation. We start with quantizing the higher dimen- 
sional field theory in the standard way of QFT. We require the dominant configu- 
ration, the path r = r\y(y), to be equal to the solution of the minimal area line 
r g (y) by introducing the weight function W(p, y). The bulk geometry takes part in 
the "coordinates quantization" by fixing the central configuration in this way. The 
closed-string-like configuration comes into this formalism through the minimal area 

Ref. 24 ' where the idea of the phase space (y,p) was presented in relation to the divergence problem 
of the "deformed" propagator in the 5D bulk-boundary theory. 
*' This situation is similar to that in the lattice gauge theory. 



32 



S. Ichinose 



principle. The field theories, which are applicable to this formalism, are limited to 
those that are renormalizable on the 3-brane, such as the 5D <£ 4 -theory, 5D YM 
theory, and 5D QED. However the advantages of this approach is that it is based on 
the QFT, hence we can expect various phenomenology applications. We stress the 
practical importance of calculation. 

We have pointed out a possibility of quantizing higher dimensional field theories 
within QFT. 
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Appendix A 

Minimal surface curve in 5D flat space 

We analytically examine the minimal area equation (15-60 in the 5D flat geometry. 

3 -TT-2 =0 ' °<V< 1 > (AT) 
1 + r z 

where r = dr/dy and r = d 2 r/dy 2 . 

A.l. Classification of minimal surface curves 

Before solving the equation, we classify all solutions using the differential equa- 
tion above. This is useful in drawing minimal surface curves and confirming the 
noncrossing of curves. (The requirement comes from the renormalization-group flow 
interpretation of the curves. See a few sentences before Eq. (j5-5|) .) 

In terms of u = 1/r 2 , the above one (Eq. (|ATp ) can be expressed as 

^) = ^ = ^r>° J n = -6n 2 <0 , 0<y<l . (A-2) 
From this equation, we can derive an important inequality relation: 

u\ y= i — u\ y =o = —6 / u 2 dy < . (A-3) 
Jo 
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The inequality u < in Eq. (IA-2j ) implies that u(y) is convex upwards. 

Making use of the above relation, we can classify all solutions (paths) as follows. 

(i) u{y = 0) > 

(ia) u(l) > 

In this cae ii(y) > for < y < I. u(y) is simply increasing (r(y) is simply decreasing), 
Figm 

(ib) < 

(iba) u(0) < u(l) Fig. [J7] 
(ib/?) n(0) > u(l) Fig. UM 

(ii) u{y = 0) < 

u(y) is simply decreasing (r(y) is simply increasing), Fig. [TUl 



5 

4 . 5 

4 

3.5 

3 

2 . 5 

2 

1 . 5 

0.2 0.4 0.6 0.8 1 

Fig. 16. Geodesic curve r(y) of Eq. (|A-l[l by Runge-Kutta method. Type (ia) simply decreasing. 
r(0) = 4.472, r(0) = -22.36 




2 . 4 
2.2 

2 

1 . 8 
1 . 6 
1 . 4 
1.2 

0.2 0.4 0.6 0.8 1 

Fig. 17. Geodesic curve r(y) of Eq. l|A-l[) by Runge-Kutta method. Type (iba). r(0) = 
2.236, r(0) = -5.590. 




Although numerical solutions are displayed in Figs. [T6rll9l the general analytic 
solution is given in Appendix A. 2. We have confirmed the high-precision equality 
between the numerical curves and the analytical ones. 
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0.2 0.4 0.6 0.1 



Fig. 18. Geodesic curve r(y) of Eq. (|A-l[) by Runge-Kutta method. Type (ib/3). r(0) 
1.4142, r(0) = -1.4142. 
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Fig. 19. Geodesic curve r(y) of Eq. ()AT[) by Runge-Kutta method. Type (ii) simply increasing 
r(1.0) = 10.0, r(1.0) = 350.0. 



A. 2. General analytic solution of minimal surface curve (Eg. h5- 6\) , (HHP, or II A- Sty ) 

We solve the differential equation (Eq(refsurf2), (|A-ip . or (|A-2|) ). which is the 
minimal surface trajectory of the 5D flat geometry. The first integral is given by 

d (\ 



aly \2 



— [ -u z + 2u- 







1 



—u 2 + 2n 3 = 2C , C(> 0) : an integral constant 
da 



u = ^ = ±2VC - u 3 , 0<u<VC 
ay 



(A-4) 



We see that the present classical system is equivalent to one particle mechanics with 
the potential V = 2-u 3 . The second integral is obtained as follows: 

J —j=== = ±2 J dy = ±2y + C', C' : another integral constant , (A-5) 

The LHS of the above equation can be integrated as 
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where F(ip, k) is the elliptic integral of the first kind and the following integral formula 
is used. 



dx 



F{ip,k) 



-^F \ cos 



.! V3-1 + V \/3 + l 



u<^- 1 + °<l , 



V3 + 1 



d0 



Vl-fc 2 sin 2 ./o x/(l-z2)(l-^ 2) 



= F( Zl ,k) , 



<k<l , _ IT < < f , 01 = sin<£ = sn(F, fe)(A;7) 

where sn(.F, k) is defined by the inverse function of F(zi,k) and is called Jacobi's 
elliptic function. It has the period AK{k) = 4F(^,k). Hence, 

Z\ = sinip = ±\/ 1 — cos 2 if = ±{(1 — cos</j)(l + cosy?)} 1 / 2 



V3 + 1 



(A-8) 



We can solve the above equation w.r.t. -u = ■^=. E3 



1 3/^(v / 3 + ^sn 2 .- 2v/3(l - |cn+|) 
U+ = o = VO ^ 

r+ 2 sn 



+ 



1 3/ -(v / 3 + l)sn 2 - 2^(1 



cm 



sn_ 



m±(y,C,C') = S n (-VsC^^y + C"), ^j=f\ 
cn ± (y, C, C) ee cn (-^C^(±2y + C"), ^-^j 



(A- 



From the requirement u = \ > 0, the above result indicates that I**) 



y/3-1 
\/3 + l 



<|cn ± (y,C,C")| 



cn -^3C 1/6 (±2y + C7' 



„ \Z3 + l ' 
2^ 



< 1. (A- 10) 



We must choose C and C in such a way that the above relation is valid for v y E (0, 1). 
r(y = 0) and ^ are related to C and C" as follows: 



y=0 

r(0) = C~ 1/6 



|snO(C, C") 



{(V3 + l)snO(C, C") 2 - 2\/3(l - |cnO(C, C')|)} 



1/2 ' 



*-* Note that one of two solutions (of a quadratic equation w.r.t. 11) is omitted here because of 
(V3 + l)sn 2 -2^/3(1 + |cn|) = — {(>/3 + l)|ch| + l}{|ch| + 1} < 0. 

Using the relation sn 2 = 1 — cn 2 , we know that (\/3 + l)sn 2 — 2\/3(l — |ch|) = {y/3 — 1 — 
(V^+I)|cn|}{|cn|-1}. 
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snO(C, C) = sn (-^C^C, ^J- 
0(C,C) = cn (-nC^C\^±J^ 



cn 



dr 
dy 



Tr(0)\C- ' 



,=o W V r(0) 6 ' 

r + (0 = <rv , 

{(v/3 + l)sn + (Z, C, C') 2 - 2^3(1 - |cn + (Z, C, C>)\)} 1/2 

r_(/) = C- 1 /6 |s " n - ( ^ C ' C/)l — .(A-ll) 

{(V3 + l)sn_(Z,C, C') 2 - 2^3(1 - |cn_ (Z,C, C')|)} 7 

The two integration constants come from the second-derivative differential equation 
(|A-ip . In the numerical approach of the Runge-Kutta method, the constants are 

given by r(y = 0) and & . It is important to check, in the simple case of the 

flat model, that the numerical solution correctly produces the analytic one. In the 
warped case, we rely only on the numerical method. 

Appendix B 

Weight function and minimal surface curve 



In § [71 we have presented a specification of the weight function W(p,y). Here, 
we examine its validity by numerically evaluating the consistency equation, to be 
given later, for each W that appeared in Eq. (|6-ip . 



B.l. Reciprocal space 

The 4D boundary manifold described in § [5] is characterized by the metric, 



ds 2 = [S ab + — dx a dx b = g ab {x)dx a dx h , r = Vx^ , 
(5 ab ) = diag(l,l,l,l),A = ( v / detg ab d A x = I Vr' 2 + 1 r 3 dy, r' = f ,(B-1) 

J Ji/a d v 

where {x a ;a = 1,2,3,4} are the coordinates of the 4D Euclidean space manifold. 
We introduce the reciprocal coordinates {p a } defined by 

x a = ^ , p 2 = p a p a ; p a = ^ , x 2 = x a x a ; x 2 = \ . (B-2) 

pA X Z pZ 

The metric (|BT|) can be, in terms of these coordinates, rewritten as 

dx a = \(-5 ab + 2Q ab )dp b , Q ab = S ab -^- , 

pA pZ 

Q ab = Q ba , p a Q ab = o, Q ab p b = 0, Q ab Q bc = Q ac , 5 ab Q ab = 3, 



ds 2 
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J^(^ab+P 2 ^jdp a dp b = g ab (p)dp a dp b , p=^¥v^ ■ (B-3) 



The 4D volume (the area of the boundary surface) is given by 



l/yl V VP / d y 

where the S 3 property: p = l/Vx 2 = r(y)~ 1 = p{y) is used. A is the UV cutoff for 
the 4D momentum integral. The minimal area principle gives the equation, 

p(y) -> P(y) + $p(y) , SA = o , 

, (p"p-2p>*)/p* d 2 p IfdpV 3 

3+ i + (pVp 2 ) 2 ° r ^ + ^UJ + P = • ( } 

This equation is the same as Eq. (|5-6p expressed by r(= 1/p). 

B.2. Numerical confirmation of the relation between weight function and minimal 
surface curve 

The dominant configuration (path, p = p(y)) for Casimir energy Eq. (l6-ip is given 
by an ordinary variation method (minimal "action" principle). 



.3 n n/. n \i c f dp d d \ 3 



y^y + 5y , 

5{W(p(y),y)p 6 (y)F(p(y),y)} = 5y (jg^ + (p 3 W(p, y)F(p, y)) = 0, 

dln(WF) 

^ = (B-6) 

dy 3 , d\n(WF) ' ^ DU ^ 
p <9p 

where W(p, y) is the weight function 'practically' introduced in §[6] and is properly 
defined in §[7] using 5D coordinates (r, y). In this subsection, we define W(p, y) using 
(p, y). F(p, y) is given in Eq. (|4-2p . The coordinate version of the above result is given 
in the text (Eq. ([7-2p ). We can also obtain d 2 p/dy 2 . 



W y = d y ln{WF) , Wp = dp \n(WF) 
Vi(W,dpW,dyW;p,y) 



dp _ Wy 



dy + Wp 

d 2 p OyWy Wy(dpWy + OyWp) ^ ( " 3p~ 2 + OpWp) 



+ 



dy 2 3p-! + Wp (3p-! + Wp) 2 (3p-! + Wp) 3 

= V2(W,dpW,d y W,djW,d y dpW,d 2 W;p,y) . (B-7) 

Note that the RHSs of the expressions dp/dy and d 2 p/dy 2 above are functions of 
p and y. If we consider that W(p, y) is unknown, by inserting the above ones in 
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Eq. (IB-5p . we obtain a partial differential equation for W(p, y) involving up to 2nd 
derivative. We call it the "W-defining equation" . 

W-defining Equation : V 2 + -V? + 3p 3 = , (B-8) 

P 

where Vi and V2 are defined in Eq. (|B-7j) . We consider that the 2nd-derivative differ- 
ential equation defines the weight function W(p,y). It is difficult to solve it. Here, 
we are content with a numerical consistency check. 

If we take some example of W(p, y) appearing in Eq, (|6-ip . the dominant con- 
figuration pw = Pw{y) is graphically shown by the valley-bottom line of Casimir 
energy integrand W(p, y)p 3 F(p, y), (see Figs. [T2HT5]) . On the other hand, there ex- 
ists another path p g (y) that is the minimal surface curve of the bulk geometry, that 
is, Eq. (|5-6p or (|B-5p . p g (y) is determined by the 5D metric, and is completely in- 
dependent of both the weight W and the model F. The VF-defining equation 
above represents the equality pw(u) = Pg(y)- We can numerically compare Pw{v) 
and p g {y). We show, in Fig. [201 the contour of Fig. [15] (parabolic suppression 
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0.1 0.2 0.3 0.4 0.5 



Fig. 20. Contour of p 3 We(p, y)F(p, y) (parabolic suppression 2, Fig. [T5)l . A = 10, I = 0.5. Hori- 
zontal axis: 1.001/A < y < 0.99999Z; vertical axis: l/l < p < A. 

W6). We can see the valley-bottom line as pw{y) ~ 4.3^. In Fig. EH we show the 
minimal surface curve p g {y) = l/r g (y), which is the r_(y) solution with C = 5.1215 
and C = 1.068 in Eq. ([A-9p of Appendix A. 2. The two graphs are similar, at least 
in shape and magnitude order. For other weights, we confirm a similar situation. 



*' Precisely, the boundary conditions at y = and y — I (boundary values: p(0) and p(l)) are 
also necessary. 
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Fig. 21. Minimal surface curve l/r_(y),C = 5.1215, C = 1.068 in Eq. (|A-9[) . Horizontal axis: 
< 2/ — 0-5; vertical axis: < l/r_ < 3. 



Appendix C 

Numerical evaluation of scaling laws: Ec as &4'5\) > ^Cas U£M> an d ^Cas 



In the text, (regularized) Casimir energy is numerically calculated in three ways: 
1) original version (rectangle-region integral), 2) restricted-region integral ( Randall- 
Sundrum type ), 3) weighted version. The final expressions show the scaling behav- 
iors about the boundary (extra-space) parameter I and the 4D momentum cutoff A. 
The results are crucial for the present conclusion. Hence, we expalin here how the 
numerical results are obtained. 

First, let us take the unweighted case with the rectangle integral region (original 
form) of Casimir energy (|4- 2[) . 

2\ 2 E Cas (A,l)= [ A dp f dy p 3 F(p,y) , (C-1) 
Jl/l Jl/A 

where p 3 F(p,y) is explicitly given in Eq. (|4-4p . The integral region is graphically 
shown in Fig. [T] as a rectangle (e = 1/A,fx = l/l)- The graphs of the integrand of 
(gH|l, p 3 F{p,y), are shown for (I, A) = (1,10) [Fig. [2], (1, 100) [Fig. [3], (1, 1000) [Fig.gj, 
in the text. From the behaviors we can expect that Ecas(A,l), ()C-ip . leadingly be- 
haves as I A 5 , because the depth of the valleys, shown in Figs. [2HH is proportional to 
A 4 and the graph-behaviors are monotonic along the y-axis (except near the bound- 
aries y = 1/A and I). It is confirmed by directly evaluating Eq. (|C-ip numerically 
(the numerical integral in Ref.31)). We plot the numerical results for various A and 
I values in Fig. [221 From the straight-line behavior we can safely fit the curve as 
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Fig. 22. Casimir Energy -Ecas of (|CT[) for various (A, Z) values. Horizontal axis: ln/1 (A = 
10,20,40, • ■ • , 1280); vertical axis: — ln(|2 3 7r 2 _Ecas|)- The results are grouped into three lines. 



The values placed on the top, middle, 
respectively. 



and bottom lines correspond to I = 10, 20, and 40 



2 3 ir 2 Ecas = lA 5 (az, + b^ln^A)). The best fit is given by (manipulating numerical 
data in Ref.31)) 

2 3 -K 2 E Cas (A,l) = -0.1249 M 5 - (1.41,0.706,0.353) x 10~ 5 I A 5 ln(M) . (C-2) 

The triplet results correspond to I = 10, 20 and 40. The first term is firmly fixed 
(the number of significant figures (NSF) is 4), whereas the second term is unsta- 
ble (the coefficients are proportional to 1//). The second one is numerically small 
compared with the first, and its determination requires more careful treatment of 
small numbers. To determine it firmly, calculation using larger values of A and I is 
necessary. 

In the restricted region case, Eq. (|5Tp . we do the numerical integral of the fol- 
lowing expression, 

/•A /■!/<? rl r^/ u 

2\ 2 E™ S = dq dyq 3 F(q,y)= du dpp 3 F(p,u) . (C-3) 
Jl/l Jl/A Jl/A Jl/l 

We plot the results in Fig. [23] for various values of A and I. From the straight line 
behavior, we can safely fit the curve as 2 3 t: 2 E^ s = A i (a^ + b±hi(lA)). The best fit 
is given by 

2 3 7T 2 E§ s as (AJ) = 

-8.93814 x 10~ 2 A 4 + (+7.73, -4.83, +5.00) x 10~ 10 A 4 ln(M) . (C-4) 
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Casimir Energy Eq^ s of (|C-3|) for various values (A, I). 



Horizontal axis: ln/1 (A 



Fig. 23. 

10, 20, 40, • ■ • , 1280); vertical axis: - ln(|2 3 7r 2 Scf 3 |). The results are placed on a straight line 
for different A and overlap (within the presented dots) for three different values of I = 10, 20, 
and 40. 



The triplet results corresponds to I = 10, 20, and 40. The first term is firmly fixed 
(the NSF is 6), whereas the second term is unstable. The second one is numerically 
very small compared with the first, and we may say that the second term vanishes 
within the present numerical precision. 

Finally, we explain the weighted case in Eq. (|6-ip taking the elliptic type W\ as 
an example. 



A A 

dp / dy fW\ (p, y)F(p, y) 
l/l Jl/A 

1.557 

Ni 



1 



iVi 



8vr 2 



(C-5) 



where the UV cutoff A and IR cutoff I are introduced to see the scaling behavior. 
In Fig. [241 we show the numerical results of E^ s (A,l) x Ni for different A and I 
values. (Note that the axes of Fig. [25] are on a linear scale, not on a log scale.) From 
the straight line behavior of Fig. |2"4"[ we can safely fit the curve as 2 3 7r 2 .E' ( K^, xJVi = 
(A/l)(ai + bi hx{lA)). The best fit is given by 



2\ 2 E^ as {A,l)xN l 



•(3.892,3.894,3.894) 4 + 
l 6 



-0.221,1.70,1.76) x 10" 



A In 



The triplet results correspond to I = 10, 20, and 40. The first term is firmly fixed 
(the NSF is 3 at least), whereas the second term is unstable. This data shows that 
the coefficients become stable as I increase. As for other types of W, the best fit 
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Fig. 24. Casimir Energy E^ s of Q&ty for various [A, I) values. Horizontal axis: A (A — 
10, 20, 40, • • • , 1280); vertical axis: 2 3 7r 2 

E cL x ^i- The results are grouped into three lines. 
The values placed on the top, middle, and bottom lines correspond to I — 40, 20, and 10 respec- 
tively. 



scaling behaviors are listed in Eq. (|6-4p of the text. The behaviors taking W2 and W3 
are different from the others. This is because the normalization factors N2 and N3 
are divergent. The leading values, except W2 and W$, do not so much depend on 
the choice of W. 
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